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Abstract 



■ A procedure for the evaluation of correlators of any order in a rea- 

^ , sonable computer time is presented. Connection between correlators 

and fluctuations of the event mean values of observables is discussed. 
Extension of the procedure to event-by-event approach is suggested. 
The usefulness of the method is demonstrated using the events simu- 
' lated within various models of multipaticle production. 

o 

°. 

Keywords: particle correlations, correlators, mean value fluctuations 

in 
o 

1 Introduction 

o 

Recent and near-future experiments devoted to the study of relativistic nu- 
clear collisions (see, e.g., [1, 2]) have stirred up considerable interest in suit- 
able methods of the analysis of high particle multiplicity events. For example, 
$_i ■ fluctuations of various observables are widely studied using the event-by- 



event approach [3]. In this paper we consider the integral characteristics of 
particle correlations, namely the correlators [4] for a given particle observ- 
able (e.g., the particle energy E, transverse momentum p t or rapidity y). 
We suggest a fast procedure how to construct these quantities. We also in- 
vestigate their properties, including the connection with the event-to-event 
fluctuations, and their ability to reflect underlying mechanisms in the case 
of several models. 

The paper is now organized as follows. In section 2 we give the formal 
definition of correlators and discuss some of their properties. Sections 3-5 
deal with the evaluation of correlators from experimental data. In section 4 
we also consider the connection of the correlators with the fluctuations of the 
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observable event-mean values. In section 6 we apply the proposed procedure 
to events simulated within various models of particle production. The results 
are summarized in section 6. In Appendix we derive the expressions allowing 
one to relate the correlators with the central moments of the single-particle 
distributions and the event-to-event fluctuations of the observable mean. 

2 Particle correlators 

Let us start with events of fixed multiplicity n of the produced or observed 
particles of a given kind. The inclusive or exclusive production probability 
of v < n identical particles can be described by a multivariate probability 
distribution function (PDF) fj) n \xi, . . . ,x u ) of a given observable x, where 
ff) n ^ is a symmetric function normalized to unity: 

/ f!) n \ x i: ■ ■ ■ ,x v )dxi . . .dx v — 1. (1) 
The PDF can be characterized by the mean value 

x {n) = J xf[ n \x)dx, (2) 

and - by the /-order moment-type quantities 

Cfci = / i x n ~ * (n) ) ■ ■ ■ K - x (n) )li n \xu • • • , x n )d Xl . . . dx n . (3) 

Due to the PDF symmetry, all the quantities C}"^ coincide for any permu- 
tation of indexes. 

In the case of independent production of particles of a given kind the 
multivariate PDF reduces to the product of the same one-particle distribution 
functions (neglecting the correlations due to quantum statistics and final 
state interaction): 

ft\x l ,...,x u ) = jf\x l )...ff\x u ). (4) 

Obviously, in this case C^fJ A = for %\ ^ %i . . . ^ %\. The quantities 

C^fJ ^ for ii 7^ . . . 7^ %i thus measure a correlation among produced 
particles and therefore we will call them correlators. Since correlators do not 
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depend on the particular set of I mutually different indexes, we will use for 
them the simple notation c[ n,x \ Note that for equal indexes i\ — %2 — ■ ■ ■ — 
ii, C^fJ [ A = M^'^ represents the /-th central moment of the one-particle 
PDF.' 

In practice, number n of produced or observed particles of a given kind 
can vary from to the ultimate multiplicity n max , so one has to introduce 
the corresponding probabilities p n and rewrite Eqs. (2) and (3) as 

flmax ^max ^max p ^max 

(x) = (iW) = E PnX (n) / E Pn = E Pn / xf[ n \x)dx/ E Pn, (5) 
n=l n=l n=l n=l 

^max 

= E Pn' 

n=l 

/^max 
(xi - (x)) . . .(xi - {x))ff\x u . . .,xi)dx x . ..dxi/ E Pn- (6) 

n=l 

Generally, the PDF can depend on various event characteristics a, includ- 
ing the observed multiplicity, the particle composition or the selected range 
of the impact parameters. Eqs. (5) and (6) should be then generalized by 
the substitution J^Pn — > J2 Pa- Particularly, 

(x) = (xM) = 5>aX< a >/5>« = E^ / xfi a \x)dx/J2p a . (7) 

a a a a 

Using the identity 

x k - (x) = (x k - x™) + (x (a) - (x)), (8) 

one can rewrite the global correlator in terms of the a-dependent cor- 
relators 

C[ a ' x) = J( Xl - (*<">» . . . (x, - (xW»/, (a) (xi, • • • , Xl )d Xl ...dxt (9) 

and the fluctuations of the observable mean x^ at a given a around the 
global mean (x): 

c\ x) = Y,p* E ([)c ( r\x (a) - (^))'- A /E^ 

a A=0 V / a 
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^(eQ^^-^), (io) 

where C^ x) = l,c{ a,a ° = 0. One may see that the absence of the correlation 
at any a (i.e. cf a '^ = for I > 1) does not lead to a vanishing global 
correlator. In this case, the latter is solely determined by the fluctuation of 
the a-dependent observable mean: = ((x^ — (x)) 1 ). 

The above formalism ignores the possible non-identity of the selected 
particles. Its generalization to the correlators of different particle species is 
however straightforward. For example, for two types of particles, say those 
characterized by positive (+) and negative (— ) charge, one has to make the 
substitutions I — > l + , /_, x — > x + and x_. Particularly, 

x\^ = J x + f[J(x + )dx + , x^ = J X-fQi{x-)dx-, (11) 

ct + t = J (s+i - (4 a) » • • • (*+<+ - (4 a) »(^-i - $- a) » • • • (*-,_ - <^ a) »- 

fi + l-( x +n ■ ■ ■ i x +i + 'i x -n ■ ■ ■ ■> ^- L )rf^+i • • • dx +l+ dx_ 1 . . . dx_ t _ . (12) 

It is instructive to express the correlator = of two charged 
particles in the events with n = n + + n_ selected particles through the 
correlators eg = Cg, eg = C^l and = Cf2. Using the identity 

x = x± — (x± — x) = x± — Ax±, (13) 

one can write 

c g = !^ [c W +A ^ + 2 ] + !h^[ C W+Ax_ 2 ] + ^[c2 + Ax + Ax_], (14) 

T^cc ^cc ^cc 

where 

n cc = n ++ + n + n + _ = n + {n + — l)/2 + n_(n_ — l)/2 + (15) 

is the number of charged pairs. One may see that even in the absence of 
correlations of particles of given charges, i.e. C+l = C^l = C+2 = 0, 
the correlator can be non-zero provided x + ^ X-. In particular, for 

n + = n_ = n/2, 



and, in the absence of correlations, 



<*' = J ^0W 2 (17) 

3 Correlator estimates 

Let us first consider n ev t experimental events with a fixed multiplicity n. The 
mean value of the observable x can then be estimated as 

1 ™cvt 

<*> = — £*<•■>, (18) 

"cvt j =1 

where is the estimate of the observable mean in the i-th event: 



x 



(0 _ I ^0 



E4- (is) 



Similarly, the correlator can be estimated as 



i n-cvt 

Cf^-E^- (20) 



^evt 



=1 



where C^' x ^ is the estimate of the correlator in the i-th event: 

Q (i,x) = ^E(^-<*»-(^-<*»; (2i) 

the sum in Eq. (21) runs over all n\ sets of / particles chosen from n particles 
in the event. One can sum either over ni = (j^j = n(n — 1) . . . (n — I + 
of the sorted sets i\ < i 2 < ■ ■ ■ < %i or over ni = n(n — 1) . . . {n — I + 1) of the 
unsorted sets i± ^ i 2 ^ ■ ■ ■ ^ i\- In the latter case, each of the unsorted sets 
gives rise to the /! of identical terms in the sum. 

In the case of a mixture of events with different multiplicities one can 
estimate the observable mean and the correlators in a similar way as in the 
two-step averaging procedure given in Eqs. (18)-(21). One should only make 
substitutions n — > n^' l \ rti — > nf 1 in Eqs. (19), (21) and take into account 
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that the single-event averages enter into Eqs. (18) and (20) multiplied by the 
weights proportional to and n\ l \ respectively: 

1 Wevt "evt 

(x) = (x®) = ^7 £ w(<) * (0 > ^ = E ( 22 ) 
^ i=i i=i 



j n ev t ricvt 



Cf° = (C^), = ^ E ^ = E (23) 

v ' i=i i=i 



The same result can be obtained by averaging simply over all iV collected 
particles or all iVj /-particle sets (formed from particles within the same event 
only): 

1 w 

<*> = tfE** (24) 

C^ = i-E(^-(^---K-^))' (25) 

where the sum in Eq. (25) runs over all iVj sets. 

The generalization of Eqs. (18)-(25) to two or more particle species is 
straightforward. For example, for two types of particles characterized by 
positive and negative charge, one has to make substitutions I — > l + , n — > 
n + ,n_, ni — > iVj — > an d ^ — ► Particularly, 

= (x? - <*+»(^ } - (^-)), C%f_ = C^ct-K (26) 

To calculate the errors, we can split all the events into n g subgroups, 
each with about the same number of events, and estimate the correlator 
(we omit all other indexes for the sake of simplicity) in the m-th subgroup 
using the global estimate of the mean value (x). The global estimate of the 
correlator is then given by the mean of the group values: 

1 n s 

C = — E^ M - (27) 

U 9 m=l 



The dispersion of the group values is 



D = 



n g 



1 £ (cM _ C7)2j (2g) 

1 m=l 



6 



where the sum of the deviations squared is divided by n g — 1 since one degree 
of freedom is used to determine the global estimate of the correlator according 
to Eq. (27). Neglecting a small correlation of the group correlators due to 
the use of the global estimate of (x), one can calculate the error in the global 
correlator estimate as 



1/2 

a(C) = (D/n g f 2 = ( 1 E (C {m) -Cf) . (29) 

n g \n g i) m=1 j 



4 Correlators and event-to-event fluctuations 

To relate the correlator with the event-to-event fluctuation of the ob- 
servable single-event mean one can use the analog of the identity in 
Eq. (8) with the substitutions Xk — > x^ and x^ — > x^\ and rewrite the 
estimate of the single-event correlator in Eq. (21) in the form 



cr = j:( x )cr\x® - (x)y-\ (so) 



A=0 

where Cq' x ^ = 1, = and cf'^ for I > 2 is defined similar to Eq. (21) 
except for the substitution (x) — > x®: 



c 



^^^tE^?-^)-^-^)- ( 31 ) 



1 „(0 



n,; 



The estimate of the correlator can then be written in the form 

^^(EQ^^-^r) , (32) 

where the /-dependent averaging is defined in Eq. (23). Introducing the 
notation 

Ax (i) = x {t) - (x) (33) 
and omitting the event indexes, Eq. (32) particularly yields: 

C { 2 X) = (c? + Ax') 2 , (34) 

= (d x) + Sc^Ax + Ax 3 ) 3 , (35) 
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Cf ) = ( c f + 4c? Ax + 6c? Ax 2 + Ax 4 > 4 . (36) 
For two different particle species + and — , Eqs. (26) and (32) yield, e.g.: 

= Cfl = (Ax + Ax_) 1;1 , = CS- = ((cr +) +Ax + 2 )Ax_) 2jl . (37) 

The meaning of quantities c? is clarified in Appendix, where it is shown 
that they can be expressed through the estimates of the moments of the 
single-particle ^-distribution in a given event: 

1 n 

mf = -Y.{x - x)\ (38) 

71 3=1 



Particularly, 



c? = -^r, (39) 
n — 1 



cf> = ^ (40) 
(n — l)(n — 2) 



» 



3// (m?) — 6m? 



(n-l)(n-2)(ra-3)' 
One may conclude from Eqs. (34)-(37) and Eqs. (39)-(41): 



(41) 



Quantities c? are determined by the shape of single-particle ^-distribution 
and by multiplicity of detected (n) or selected (n — > v < n) particles 
in a given event. Therefore, they are sensitive only to a part of the 
correlation related to this shape; an example is the correlation due to 
energy-momentum conservation (see the 4-th item). The remaining 
part is contained in the event-to-event fluctuations of the observable 
mean in accordance with Eqs. (32)-(37). 

The magnitude of the quantities c?^ decreases with the increasing 
number v of selected particles. This decrease should be compensated 
by the //-dependence of the fluctuations of single-event mean observable 
to guarantee the //-independence of the correlators. Particularly, for 
uncorrelated production (C^ = 0), one gets (Ax 2 )? = (m?) 2 / 'iy — 
!)■ 
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• Assuming the moments mf of the single-particle distribution weakly 
varying with the number n of observed particles, the correlators for 
high enough n are dominated by the event-to-event fluctuations, i.e. 

• For a conserved additive observable x (e.g., particle energy), consid- 
ering particles of any kind and assuming that all n tot particles are 
observed (n = n tot ), the mean x does not fluctuate so = (cf^)j. 
Particularly, Cf } = -(m^) 2 /(n tot - 1) and C { 3 x) = 2(m^) 3 /[(n tot - 
l)(^tot — 2)] for the events with about the same numbers of produced 
particles. 

• Evaluation of correlators c\ x ^ using Eqs. (39)-(41) can save substantial 
amount of computing time compared with the direct evaluation accord- 
ing to Eq. (25). In the former case the number of operations is oc n 
while in the latter case it is can 1 . 



5 Event-by-event correlators 

In the case of uncorrelated particle production, the event-to-event fluctua- 
tions of the observable mean are solely determined by the quantities cf\ 
Such "standard" fluctuations are thus related to the single-particle distribu- 
tions and can be estimated by the event-mixing techniques. The eventual 
deviation from the "standard" or "statistical" fluctuations then signals the 
presence of correlations or - a mixture of the events with different single- 
particle distributions. 

Therefore, to clarify the origin of the "non-standard" or "non-statistical" 
fluctuations, it is desirable to estimate the correlators for the events with sim- 
ilar characteristics a. The ultimate solution is to use the information from 
a given event only. To do so, one has to destroy the equality in Eq. (A.l), 
i.e. decouple as much as possible the observables Xj entering into the corre- 
lator estimate from the observable mean x in a given event. In the case of 
sufficiently large multiplicity, it can be achieved by splitting the event in a 
number of sub-events s = 1,2, .. . ,n scv t with about the same multiplicities 
n^ s K One can then estimate the correlators according to Eqs. (18)-(25) or 
(34)-(37) and (38)-(41), making the substitutions i — > s and n evt — > n sevt . 
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Table 1: The correlators corresponding to PDF's in Eqs. (42) (the first two 
lines) and (43). The results of numerical integration are compared with the 
estimates obtained from n evt simulated events. The allowed interval of the 
variable x was 0-150, the parameters T = 150, Ti = 25 and A = 0, 1, 2. 



Correlator 


Cf 






n = 3 


1312.5 


0. 




n ev t = 5 • 10 5 


1313.4 ± 1.8 


61. ± 108. 




n = 3, A = 


0. 


0. 




n evt = 5 • 10 5 


-0.3 ± 1.4 


21. ±70. 




n = 3, A = 2 


229.0 


1771.5 




n evt = 5 • 10 5 


229.9 ± 1.4 


1779. ± 110. 




n = 4, A = 1 


113.8 


707. 


2913. 


n evt = 3 • 10 6 


114.4 ±0.5 


666. ± 23. 


3890. ± 1836. 



6 Examples 

Let us first consider the following simple PDF's characterized by the param- 
eters T,T\ and A: 

/r^i^^oce-'^'e-' n l e "l l) ? (42) 

ti n \ Xl ,x 2 ,...,x n )cxf[e-%(l + A ^ e" 1 ^ 1 ), (43) 

The parameter A in Eq. (43) controls the correlation strength; for A = 
there are no correlations among particles, so = 0. Eq. (42) corresponds 
to two-particle correlations only, thus leading to = 0. 

The correlator expectation values calculated according to PDF's in Eqs. 
(42) and (43) by direct numerical integration and their estimates obtained 
from simulated events are compared in table 1. One may see that the esti- 
mated values of the correlators agree with their expectation values within the 
errors. The relative errors rapidly increase with the order / of the correlator 
thus making quite uneasy its measurement for / > 4. 
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Table 2: Correlators C\ and C\ estimated from the "microcanonical" 
ensemble of events with fixed mean particle energy x = 100 in arbitrary 
units, each consisting of 2 • 10 4 events with fixed particle multiplicity n tot . 
All particles are assumed to be observed, i.e. n = n tot . 



n tot 


5 


10 


15 


20 


100 




-1181 ±2 


-624 ± 1 


-426 ± 1 


-323 ± 1 


-66.2 ±0.1 




49900 ± 200 


14730 ± 50 


6940 ± 20 


4041 ± 6 


174 ±1 



Table 3: The same as in table 2 for n tot = 100 and different numbers v of 
selected particles. 



V 


5 


10 


15 


20 




-75 ±8 


-63 ±3 


-68 ±2 


-66 ±2 




200 ± 500 


300 ± 100 


220 ± 70 


170 ± 50 



As another example, we consider the energy correlators estimated from 
the "microcanonical" ensemble of events simulated according to the non- 
relativistic phase space (non-relativistic ideal gas of particles) using Metropo- 
lis algorithm to redistribute the particle energies via binary energy-conserving 
collisions. The total amount of energy distributed to n to t particles is X tot = 
ritotx, where x is the mean particle energy. We have put x = 100 in arbi- 
trary units. Assuming that all n tot particles are observed, the mean particle 
energy does not fluctuate (Ax = 0) and the correlators are then given 
by Eqs. (39)-(41) with n = n tot . Particularly, and vanish at large 
n to t as l/n tot and l/n tot 2 respectively. The results shown in table 2 confirm 
this behavior. 

In table 3, we show the same correlators as in table 2 for n tot = 100, but 
now calculated for different numbers of selected particles v < n to t- One may 
see that within the errors the correlator estimates are //-independent. This 
confirms the conclusion at the end of section 4 about the compensation of 
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Table 4: The estimates of the pion transverse momentum correlators obtained 
from 2 • 10 3 simulated central PSM events of Pb + Pb collisions at total cm. 
energy y/s = 200 AGeV. 



Species 


(GeV/c) 2 


(jipt) 
(GcV/c) 3 


(jipt) 
(GeV/c) 4 


7T + 


(44.0 ±6.6) • 10~ 6 


(-7.4 ± 12.1) • 10~ 8 


(8.7 ±3.3) • 10" 9 


7T~ 


(44.6 ±6.0) • 10~ 6 


(7.8 ± 12.5) -10- 8 


(9.5 ±4.6) • 10" 9 


7T U 


(71.2 ±7.6) • lO" 5 


(3.0 ±2.7) • 10" 7 


(34.3 ±9.5) • 10- y 


7T± 


(52.6 ±4.2) • 10- 6 


(-1.5 ±7.1) • 10~ 8 


(9.0 ±2.0) • 10" 9 


7r ± -° 


(59.6 ±3.5) • 10~ 6 


(2.1 ±5.3) • 10~ 8 


(11.2 ± 1.6) • 10~ 9 



the //-dependence of the quantities c\ ' by that of the x fluctuations. 

Finally, we have estimated the pion transverse momentum correlators 
ci Pt \ Cg^ and C| p ^ using 2- 10 3 central (impact parameter 6 = fm) events 
of Pb + Pb collisions at total cm. energy y/s = 200 AGeV simulated in 
the parton string model (PSM) [5]. In this phenomenological model, the 
soft and semihard parton collisions initiate the formation of color strings 
and the subsequent string breaking leads to the production of stable (with 
respect to strong interaction) hadrons and resonances that are forced to de- 
cay. We have divided the simulated events into 100 subgroups and selected 
v = 1900 pions of a given charge in each event. As one can see from ta- 
ble 4, the correlators are consistent with zero within the errors while 
this is not the case for the correlators and C^. The latter appear 
to be non-zero for any combination of pion charges, being somewhat higher 
for neutral pions. We have also selected it + it~ pairs only and estimated 
C^J = (61.9 ±4.4) • 10- 6 (GeV/c) 2 . Since the PSM takes into account only 
the resonances from the lowest SU(3) multiplets that do not decay into pairs 
of like-sign charged pions, the small difference between the correlators for 
pairs of like-sign and unlike-sign charged pions indicates that the non-zero 
PSM correlators are mainly of the non-resonance origin. Also, since we con- 
sider the high multiplicity events, the correlations due to energy-momentum 
conservation are of minor importance in accordance with the discussion in 
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section 4. The probable source of the non-zero PSM correlators are thus the 
semihard parton collisions that are becoming important with the increasing 
energy and are known to lead to a noticeable "non-statistical" p t fluctuation 
in Au + Au collisions at RHIC [2] . 

7 Conclusion 

We have developed a fast procedure allowing one to calculate, in a reasonable 
computer time, the particle correlators of any order and estimate their errors. 
The corresponding C++ code is available on the request at e-mail address 
amelin@sunhe.jinr.ru. We have suggested the extension of this procedure for 
the event-by-event approach as well. We have shown a close relation between 
the correlators and fluctuations of the observable event-mean values. We 
have applied the proposed procedure to the events simulated within various 
models and demonstrated the usefulness of the two-, three- and four-particle 
correlators; the measurement of the higher-order correlators is rather difficult 
as it requires very high statistics. 
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Appendix A 

To clarify the meaning of quantities cf \ let us follow Ref. [6] and consider 
the powers of the equality 

= ±( Xj -x) = ±A r (A.l) 

3=1 3=1 
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Thus, the second, third and fourth powers of Eq. (A.l) yield 

o = E A A^ + E A /> (a.2) 



Using further the relations 



= -(ev) 2 +2ev, 



and 



(A.3) 



= £ A, A fe A, + (fj £ AjAk 2 + £ A/, 
= E A j A fe A;A m + (f) E A i A fc A / 2 + 

(A ' 4) 



S4,V=-EV H , (A.5) 



E A^A, 2 = E A,A fe E A, 2 - 2 E A ,V 

#k?4i \o+k j l j+k 



(A.6) 



EVA fc 2 =fEvV-EV, ( A - 7 ) 

one can express the multiple sums in Eq. (31) through the sums of the powers 
of Aj related to the estimates of the central moments in Eq. (38). 
Rewriting Eq. (31) in the form 



,(*,*) L 



n(n - 1) . . . (n - Z + 1) 



E A n A, 2 ...A„ (A.8) 



(x) 



one then proves Eqs. (39)-(41) relating quantities c\ with central moments 



(x) 
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